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THE ROLE OF INSTABILITY WAVES IN PREDICTING JET NOISE 


M.E. Goldstein* and L.M. Handler 1 
National Aeronautics and Space Administration 
Glenn Research Center 
Cleveland, Ohio 44135 


Abstract 

Debate over whether linear instability waves play a role 
in the prediction of jet noise has been going on for 
many years. Parallel mean flow models, such as the one 
proposed by Lilley, 2 usually neglect these waves 
because they cause the solution to become infinite. The 
present paper solves the true non-parallel acoustic 
equations for a two-dimensional shear layer by using a 
vector Greens function and assuming small mean flow 
spread rate. The results show that linear instability 
waves must be accounted for in order to construct a 
proper causal solution to the problem. 

Introduction 

Current industrial noise prediction methods, such as the 
MGB approach, 1 are based on Lilley’ s 2 equation which 
involves a wave operator appropriate to sound 
propagation on a parallel mean flow. This equation 
possesses homogeneous solutions corresponding to 
spatially growing instability waves on that flow and its 
complete solution consists of a particular solution plus 
these homogeneous contributions. But the result is 
meaningless because the instability waves become 
unbounded (infinite) far downstream in the flow. The 
usual resolution to this dilemma is to completely 
neglect the contribution of the instability waves. 
Dowling, et al. 3 attempted to justify this by invoking 
causality arguments. 

A better approach is to begin with an equation 
appropriate to sound propagation on a non-parallel 
flow, say the actual mean flow in the jet. The most 
important difference between this approach and Lilley’ s 
parallel flow result is that the homogeneous solutions to 
the acoustic equations correspond to instability waves 
that grow and then decay on the diverging non-parallel 
base flow and therefore always remain bounded, which 
eliminates the dilemma alluded to above. 

A vector Greens function approach is used in reference 
4 to obtain a general solution to the non-parallel 
acoustic equations by assuming that the spread rate of 
the mean flow is small. The relevant Green’s function, 
which must satisfy an appropriate causality condition, 
consists of two components — both of which act on the 
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same source term. The first corresponds to the usual 
Lilley equation solution but with slightly modified 
source terms. The second is associated with linear 
instability waves but is very different from 
conventional instability models since these waves are 
now continuously generated along the length of the jet 
and do not constitute separate sound sources. They 
merely serve as conduits to carry the sound to the far 
field. 

The Greens Function can be thought of as a filter that 
only allows certain parts of the source to radiate. Each 
of its two components responds to different portions of 
the source spectrum causing the radiation field to 
exhibit a bi-model structure. The contribution of the 
second Greens function component is fairly small at 
subsonic Mach numbers but can be significant at 
supersonic speeds. The result is therefore somewhat 
different from the conventional parallel flow models. 

The present paper applies the general result of reference 
4 to a two-dimensional shear layer with fully developed 
similarity form, which greatly simplifies the final 
equations. 

The Fundamental Equations 

We let 

h o= h + ^ v2 (!) 

denote the stagnation enthalpy, h the enthalpy, 
t the time, x = {x\p: 2 ^} the Cartesian coordinates, 
p the pressure, p the density, v = {vi,v 2 ,v 3 } the fluid 
velocity, cry the viscous stress tensor, and q : the heat 
flux vector. The flow variables are assumed to satisfy 
the ideal gas law 

P = pRT, h = c p T (2) 

with R = c p - c v being the gas constant, c p and c v the 
specific heats at constant pressure and volume and T the 
absolute temperature. 

References 5 and 19 show that the Navier-Stokes 
equations can be recast into the convective form of the 
linearized Navier-Stokes equations by dividing the 
dependent variables 

p = p + p\ p = p + p', h-h+h', v ( - = Vj + v'j 

( 3 ) 
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( 11 ) 


as well as the viscous terms ay and q h into their ‘base 
flow’ components p , p, h, v ; - , a y , and q t and into their 
‘residual' components pi, p', h', v/, a-j , and q\ and 
requiring that the former satisfy the inhomogeneous 
Navier-Stokes equations 


D n _ 

— p = 0 
Dt 


D __ dp d 

— pvj +— = 

Dt dxj 8xj 


( f ij +a ij) 


(4a) 

(4b) 


—pK-— = 

Dt 0 dt 


dp D 0 H 0 


Dt 


— (Hj 


d Xj 


-qj+ v iCTij 


) (4c) 


along with an ideal gas law equation of state, 


h = c p f — 
p R p 


(5) 


where the operator DJDt is not the usual convective 
derivative but is defined by 


£o_ fs ^ + ±_ 

Dt dt dxj 



(6) 


for any function/ 


1 ~2 
— V 
2 


(7) 


is the base flow stagnation enthalpy, and the ‘sources 
strengths’ Tjj,H 0 , and Hj, which are assumed to be 
localized, can be arbitrarily specified. 

The residual variables are governed by the convective 
form of the linearized Navier-Stokes (LNS) equations 


I j a\U y s p, v 1, 2, ... 5 (8) 

where the linear operator L^,, which is exactly the same 
as the operator that would be obtained by linearizing 
the Euler equations about the base flow (4), is written 
out in full in reference 5. The Greek indices range from 
1 to 5, the Latin indices from 1 to 3, the Einstein 
summation convention is being used, the five- 
dimensional dependent variable vector u v is given by 

{u v } = \pv' v , jo'j (9) 


where 


Po =P' + ^Y~Pv' 2 +[y-\)H 0 
The source strengths e' iv - e iv are given by 


y — 1 2 

4 = -P V i V j + ~ 5 ijP V ' + a ij ’ e /4 = ~P v 'i h o -<li +a ij v j 

( 12 ) 

% - %j Sy e /4 =Hi -TijVj (13) 


where 


ii it - 1 il 

n„ = n — v 


(14) 


Equations (8) are exact consequences of the original 
Navier-Stokes equations, but they have been rearranged 
so that their left-hand sides are the same as the 
equations that would have been obtained by linearizing 
the convective form of the Euler equations about the 
base flow (4). In other words, they are just the LNS 
equations with modified dependent variables and with 
the viscous stress perturbation replaced by the 
generalized Reynolds stress e\: - e y and the heat flux 
perturbation replaced by the generalized stagnation 
enthalpy flux -e,- 4 .In the acoustic analogy approach 
the right hand sides would be treated as known source 
terms and the linearized Euler operators Z^ v would be 
inverted to solve for the dependent variables. The fact 
that they are non-linear causes no particular difficulty 
since one of them, p' Q , reduces to the linear pressure 
fluctuation in the far field. 

The base flow equations (4) reduce to the usual Euler 
equations when the arbitrary source strengths /y , Flj , 

and H 0 and viscous terms a j and q t are set equal to 

zero. A general class of solutions to these equations, 
which conveniently provide good approximations to the 
actual mean flow fields in jets and shear layers, is the 
unidirectional transversely sheared mean flows 

Vj = SjflJ (jc 2 , X3 ) , P = constant, p = /5(v 2 , ) ( 15) 

The fifth equation (8) then decouples from the 
remaining four which become the inhomogeneous 
compressible Rayleigh equations 6 


and the source term s M is given by 

S M -*nj) +S K 4 (4 ~%)f“ for ^ =1 ’- 4 

( 10 ) 


p 


D n u / 


- +s nuj^- 


8 Po 8 1 ( ir , 

+ — = — e ij (16a) 
dx t dxj 


+ r p^ = ( r _l) 

Dt dxj V ’ 


r Se\j dU ' 
+ e ij 


\ dx J 


dx 


(16b) 


J J 
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where u- = u t / p and £> a /Dt now reduces to the usual 
convective derivative 


(23) 



It is well known (see chapter 1 of ref. 7) that the 
velocity-like variables uj can be eliminated between 
these equations (by taking the divergence of the first 
equation and the convective derivative of the second, 
subtracting the results and then using the first equation 
to eliminate the velocity fluctuation on the left-hand 
side) to obtain the Lilley 2 equation 

Lp' 0 =s ( 18 ) 


where 


£ o _ ±~ 2 

Dt dxj dxj Dt 2 


dU 8 ~ 0 d 
dxj dx\ dxj 


is the variable density Pridmore-Brown operator, 
c2 =yp / p(x 2 ,x 3 ) 


is the square of the mean flow sound speed and the 
source term s is written out in reference 5. As noted in 
the Introduction, it possesses homogeneous solutions 
corresponding to spatially growing instability waves on 
the base flow (15). The complete solution to these 
equations consists of a particular solution plus 
contributions from these homogeneous solutions, but 
the result is meaningless because the instability waves 
become unbounded far downstream in the flow causing 
the solution to become infinite there. The usual 
resolution is to completely neglect the instability wave 
contribution. 


A better resolution is to choose the base flow to be the 
actual mean flow of the jet. The over bars in the base 
flow equations then denotes the time average 


• = lim 

T — >oo 



( 21 ) 


where the dot is a place holder for p, v h p, and /;, and 

• = [p*l p) ( 22 ) 

denotes a Favre averaged quantity 8 for all variables 
except h 0 , which is defined by Eq. (7). Notice that 
Eq. (20) is completely consistent with the overall ideal 
gas law (2) when the tilde is defined in this fashion. 

The time derivatives drop out of the base flow 
equations and their source strengths are given by 


T ij = ~P v l v j 

H 0 =\t U ( 24 ) 

hi = -fyVj = -p h' Q v'i (25) 

They are now the ordinary RANS equations, which do 
not, of course, form a closed system. They are usually 
closed by assuming some sort of model relating the 
source terms to the mean flow variables v ( - , p, p , and 
h and their derivatives, such as Boussinesq model 8 ' 9 
for the Reynolds stresses and a similar model for H ■ . 

The most important difference between these results 
and the parallel flow result is that the homogeneous 
solutions to the LNS equations, which now correspond 
to instability waves growing and then decaying on the 
diverging non-parallel base flow, will always remain 
bounded. This eliminates the paradox alluded to above 
and the corresponding LNS equations can be used to 
calculate the radiated sound. The relevant solution 
consists of a particular solution that is driven by the 
sources (i.e. it satisfies causality and therefore the 
appropriate upstream boundary conditions) and a 
homogeneous solution that is driven by the upstream 
boundary conditions, say at the nozzle exit. 

The Green’s Function Solution 

The particular solution that satisfies the appropriate 
outgoing wave boundary conditions can be written in 
terms of the causal vector Greens function for these 
equations g va , which satisfies 11 

Ljuv gva ^ ^ CT 8(x-x')5(t-f) (26) 

and leads to the following formula for u v 

M l/ (x,f) = | J g VM [x,t\x',t')s /J (x',t')dx'dt' 

-CO y 

(27) 

where the symbol V denotes integration over all space. 

The derivatives can be transferred to the Greens 
function (see reference 4) to obtain 

00 

p'o = { { rj M (x,t\x',t')T M j(x',t')dx',dt' (28) 

—00 

where 
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and 


rjti s -(r - !) 


3x • 

\ ux j 


S4/j 


Mfi 

dx 'i 


§44 


T -T 


(0) 


+ sT- W +s 2 T 
+ bl JM 


( 2 ) 


+ ... 


(37) 


(r-D 2 

2 


5 jH 


' d S4t 
, dx '( 


dv[ 

dxj 


7 844 




J 


accounts for the propagation effects and 


MJ 


with 


and 


pv'^v)- pv'^v) 

v j 


v 4 = h 'o 


'M] 


a 4j 


(29) 

(30) 

(31) 

(32) 


when the bulk viscosity is zero. t,y, which accounts for 
the source structure, is just the generalized 
four-dimensional fluctuating Reynolds stress in the 
inviscid limit which is of primary interest here. 
Equation (28) therefore provides a direct linear relation 
between this quantity and the far field pressure 
fluctuations. 

Small Spread Rate Approximation 

Our interest here is in the two-dimensional shear layer 
with zero velocity in the upper stream. But the general 
solution to Eq. (26) is very complicated even for this 
simple case. However, high Reynolds number turbulent 
shear layers are nearly parallel and it is therefore 
appropriate to seek an asymptotic solution 4 ' 5 by 
expanding in the shear layer spread rate s. Since the 
mean Reynolds stresses will certainly dominate over 
the viscous stresses in this case, we assume that 
°ij=9j= °- 


Then the mean cross flow velocity v 3 and the Reynolds 
stresses 7j i3 and T 2fJ will vanish and the remaining 
mean flow variables will expand like 

vj =U(X,x 2 ) + sU (l \x,x 2 ) + ... (33) 

v 2 =sV(X,x 2 ) + £ 2 V^ (X,x 2 ) + ... (34) 

p = P{X,x 2 ) + £P (X \x,x 2 ) + ... (35) 

p = R(X,x 2 ) + £R (l \x,x 2 ) + ... (36) 


with similar expansions for h,Hj, and H a . Here 

X = £Xi denotes the slow streamwise variable. 
Substituting these into the mean flow equations and 
assuming that all Reynolds stresses vanish in the free 

streams shows that = 0, P = T^\ 

P and that the lowest order mean flow 

equations reduce to the usual two-dimensional 
compressible turbulent boundary layer equations. 

These results imply that the linear Euler operator L ^ 
will expand like 

v=v (0) +£ L J 1] +••• ( 3g ) 


and therefore that the vector Greens function g va has an 
expansion of the form 


(0) (1) 

Svu Sva ^Svtj + • • • 
where and g[^ are determined by 

V ( °W 0) =S fUT S(x-x')S{t-t > ) 
and 



*i l) =-v (1) 



(39) 


(40) 


(41) 


The operators and L^], are written out in detail in 

reference 4. The important thing is that the left hand 
sides of the first 4 components of these equations are 
formally the same as the left sides of the Rayleigh 
equations (16) and since their coefficients are 
independent of x 3 and t and depend on X\ only through 
the slow streamwise variable X, a local (causal) solution 
can be constructed by taking Fourier-Laplace 
transforms to obtain 


(r-i)*JS— J n, 


c o{ x 2) 


IC—OD —OO 


D, 


[k x U(x' 2 )-(o\ 

,j[k x{Xx - x 0+ , 3 ( X 3 -x()- f ,(l-f)] - o ^ ^ dwdk 


(42) 
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k = {k x ,k 3 } (43) 

is the wave number vector in the plane of the shear 
layer 

D' a =4r f° r ° = j = W (44) 


-(y-l) f) 

D'a- = -y, ■, i0J + U ( x 2)^T7 forcr=4 (45) 

Co( x 2) L “*lJ 


and the one-dimensional self adjoint Greens function 
G 0 ( x 2 \ x 2 ) satisfies 


L k G 0 = 


8( x 2- x ' 2 ) 


a 2 - ,22 

h T — + 1- - C ° - (47) 

1 dx 2 ( kjU-a > ) & 2 (k x U-(o) 

with 1c = \ k \ 2 is the one-dimensional reduced 
Rayleigh’s operator. 

G 0 is therefore given by 


layer is inviscidly unstable). The k\ -integration contour 
C/ C[ must therefore be deformed to lie below these poles 
(as shown in fig. 2) in order to obtain a continuous 
result , 12 ' 13 The poles correspond to the eigenvalues, say 
k x = K( (fl,k 3 ,X) of the Rayleigh operator (47). The 
integral over can then be decomposed into an 
integral over the real Ar r axis plus a contour integral that 
can be evaluated as the residue of the pole that crossed 
the real axis. The causal Greens function (43) can 
therefore be written as 


(/-■)g4j -~J j J 

-00 [h U(x' 2 )-Co\ 


, Jk^xi-x^+^xj-x ^ )-©(?-«')]- 

u ( 


G 0 (x 2 |x 2 ) dcodk 


■'1 1 1 


c o(*2)%+(*2) 

r ~i2 

kU(x'2)-o> ” A '(®, Ar3 ) 




w 0 _ (x^dadkj 


Inversion contour 


- w + (x 2 W_(x 2 ) 

G 0 = - — - — for x, > (48) 

A(2^) 


Above all singularities 
in co-plane 


L k w+ =0 

w ± ~ e ± ^ ±Xl x 2 — > ±00 


Figure 1 . — Inversion contour for Laplace transform. 



(51) 

c+, c_, and LL are defined in the obvious ways and 

A - W ( x 2 ) cT( * 2 )/[ k x U ( x 2 ) - cof (52) 

where W is the Wronskian of a>+ and (o is a constant 
independent of x 2 . The result for x 2 < x 2 will not be 
needed. 

Steady state solutions can only exist if the Laplace 
inversion contour (shown in fig. 1) can be continuously 
deformed onto the real axis. But G <0) possess (usually 
simple) poles in the upper half A' r plane that cross the 
real £ r axis during this deformation (since the shear 
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where w 0 ±, which are now linearly dependent, are the 
limits of w ± as k\ — > k and 

A' = lim — (54) 

k | — >A' dk | 

The first term in this result, which corresponds to the 
usual Lilley’s equation solutions that appear in the 
literature, remains bounded and is a uniformly valid 
approximation to the analogous term in the true non- 
parallel flow result. But the second terms grows without 
bound as xj - x[ becomes large (since Im k < 0) and 
therefore becomes invalid on the long streamwise 
length scale X - X . It can, however, be rendered 
uniformly valid within the shear layer by using the 
method of multiple scales, which amounts to replacing 

/[*•(», 



by 

" 1 X 

i — f /c(co,k2,X)dX 

a(x J 

where the slowly varying amplitude function A, which 
satisfies the initial condition 

A(X \X)^lasX^X (55) 

is determined by imposing a solvability condition on 
the next order solution. The detailed results are given in 
reference 4. 

But, this solution is still not uniformly valid at large 
transverse distances where Y = sx 2 = 0(1) , because the 
second order solution 14 (which is a generalization of the 
second order instability wave solution considered by 
Tam and Morris) now decays more slowly with 
X 2 as x 2 — » ± oo than the first order solution. But an 
“outer” solution can be constructed by using the WKBJ 
method to solve the reduced stationary (or uniformly 
moving) media wave equation that governs the flow in 
that region This result can be used to form a 
“composite” solution 17 that remains uniformly valid 
everywhere within and outside the shear layer. The 
details are given in reference 4. Our interest here is in 
the upper region and we therefore do not require 
uniformity below the shear layer. 

The final uniformly valid Greens function g 4tT can be 
written as 

00 

(r-l)gi^ = { e~ lt0 ^ t ~ t (x|x ';®) dco (56) 



X o =X-a(X 0 )Y (60) 

a = , K (61) 

i^ 2 + k i- ( y c ) 2 

with the arguments k 3 and co suppressed. The functions 
G a and A' are defined by Eqs. (48), (52), and (54). 

It is necessary to evaluate Eq. (57) in the region where 
R 2 = s 2 ,- 2 oo (62) 

with 



in order to in calculate the radiated sound at large 
distances from the shear layer. But applying the method 
of stationary phase and using Eqs. (48), (50), (51) along 
the asymptotic behavior of © and a shows that 

Mo) _ e icorlc + 0J . 

— > — — sm&'cos® 

( 2 c + 

[g Bct (*|x') + G /ct (x|x')] 

(64) 
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when the azimuthal polar angle 0 and the Mach number 
M are sufficiently large. The and 

G /cr (x|x') are defined in reference 4, cp is the 

circumferential polar angle and x = x(9,cp) denotes 
the unit vector in the x-direction. 

The Far Field Acoustic Spectrum 

The quantity of interest is the mean square acoustic 
pressure 

T 

pl =— \ P'o{ x ^)p'o{ x ’ t + t o) dt (65) 

-T 


^integration and the result would then depend only on 
the temporal correlation of the turbulence at any given 
point in the shear layer. This is probably the true in the 
cross stream direction, but T a j fi ii is expected to decay 
much more slowly in the stream-wise direction- 
especially at the higher Mach numbers of technological 
interest. However, Lighthill 15 pointed out that the 
cj | -decay should be much more rapid in a reference 
frame moving with the convection velocity, U c , of the 
turbulence and Ffowcs Williams 16 showed that this idea 
can best be implemented by introducing the moving 
frame correlation 

tin (*';£*■) s ?ojnt (*';£>) ( 69 ) 


at large distances above the shear layer, where we 
have used the fact that p' 0 — > p' in this region (see 
Eqs. (11), (23), and (24)) and T denotes a large but 
finite time interval. Substituting Eq. (28) into this result 
and changing integration variables shows that 

00 

p 2 (x,i 0 ) = | j jFjofje{ x \ x '’C’t 0 +r ) 

—oo V 

T a ini ( x'; £7 t ) dx'dgdr 

( 66 ) 


where 

Z = C-iU c T (70) 

into Eq. (66) and changing the integration variable from 

The result, which is much simpler when re-expressed 
in terms of the acoustic spectrum 

00 

U*) = ^ { e ia *°p 1 {x,t 0 )dt 0 (71) 


where 

00 

Yjalp = J rja{ x \ x '’h+to +T yi!fi{ x \ x ' + C’h) dt l 

—00 

(67) 

accounts for the acoustic propagation and mean flow 
interaction effects and 

T 

^ J T *i ( *'» t'Yuj { X ' + C,t' + T)dt' (68) 

-T 

is the density weighted fourth order 2-point time 
delayed turbulent velocity/total enthalpy correlation. 

Equation (66) provides a direct linear relation between 
the acoustic pressure and the fourth order correlations 
of the turbulent fluctuations in the shear layer. 
Unfortunately, the latter quantity is very difficult to 
determine — either experimentally or numerically — and 
it is, therefore, highly desirable to make the acoustic 
predictions as insensitive to its details as possible. It 
would be possible to make considerable progress 
toward achieving that objective if the propagation 
factor, Yjai.it Were to vary only slowly over the 
correlation length of ^ a ju( > because the former could 
then be treated as a constant relative to the 


or better yet in terms of J ( „ (x I x') 

I a>{ x ) = \ I <o( x \ x 'Y x ' ( 72 ) 

V 

the acoustic spectrum at x due to a unit volume of 
turbulence atx',is given by 

00 

2 k J |r 7CT (x|x>)r* jU (x|x' + ^ + iC/ c r;©) (73) 

—00 V 

e- l0>T Tcj^(x',Z,T) d %dT 

where capital letters denote Fourier transforms (with 
respect to time) of the corresponding lower case letters, 
and the asterisk denotes complex conjugates. 

When# and M are sufficiently large, equation (64) and 
the form ofGg^ and G/^ given in reference 4 imply 

that 

e l<orlc +oj . n 

r /cr = sin 0 cos cp 

(2 ;r) c + r 

[r S y CT (*|x') + r /yCT (x|x')] 
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(74) 


where 


Y Bja=^Bja{x\x' L )e 


.0) Q 

-l X\ cos t) 

C, 


= D jCTjU G BjU (x |x') 


(75) 


r ija =f /jCT (i|xi)e 


-® m (e,<pX) 


= D, 


'j<r/j G I/j ( x | x j_) 


0„ 


“o - 

: [k (Z)dZ -yi-sin 2 0sin 2 cp cosB (76) 

c + 

x’ + 


a 0 is determined implicitly by 


, tt> . . . _ , ff> . 

k(co, — sm0sin(p,a o ) = — cos0 (77) 

Cj. c. 


T ~ x d du ir-i) _ 5 

jaM ^ dx obc} CTl ^ 2 2<T SxJ, 


(78) 


K jMt {x';£,a>) = ^ j (*';£,r)rfr ( 81 ) 


denotes the 2-point fourth order turbulence spectrum 
relative to the moving frame and the second terms in 
the square brackets are identically zero when co is 
greater than the neutral frequency where 
1 m k[co x' 0 ) = 0 . 

This result is nearly exact. The only significant 
approximation is that the mean flow spread rate, s, is 
small. But is expected to vary much more rapidly 
with £ than F efl and Eq. (80) can therefore be 
approximated by 


7 r4 X l X ') = 


/ \ 2 

ocos^sin# 

c + r 2;r 


[r^-a(*|*') + r ^ (*!*')] 

T *Bin ( x | (■ x';co(l-M c 


+r 


cos#)) 

{&,*') ) 


and x' ± = {^2 , X3 } is the cross flow source variable. 


(82) 


Since 0* varies on the slow streamwise length scale X , 
which is much longer than the correlation length of the 
turbulence, we can account for its variation over the 
range of integration in Eq. (73) by expanding it in a 
Taylor series to show that 

©«, (x, X' + * (6 + U c r)) + ©» (i,*0 

-£k{X')(^+U c t) 


with more than enough accuracy to evaluate this result. 
Using this in Eq. (76) and inserting the result along 
with Eqs. (74) and (75) into Eq. (73) shows that 


f ■ a < 1 
eosvx&coscp 



c + r2n 

J [rv( x K 

b 

l-H 

+ 

V 


V*Bt M {i\x' + i 


+r *ie M {x \ x ' +i - 

’Kim 




where 


( 80 ) 


where 


( x ‘ v- ®) - J H aj ^(x';S,c0)dZ (83) 

v 

is the single point fourth order turbulence spectral 
tensor at x'. 

This approximation cannot remain valid when 
M c cos 6 = 1 because the first spectral function in 
Eq. (82) does not — » 0 as co —> 00 in this case, which 
causes the integral of Eq. (82) with respect to co (and 
therefore the mean square pressure) to become infinite. 
Ffowcs Williams 16 argued that this type of difficulty 
could be overcome by replacing the Doppler factor 
( 1 -M c cos 6) with 

\-M c cos# — > cos#) 2 +(aM c ) 2 (84) 

where a is a small constant. While this result is not 
strictly valid in the present context, it may be useful for 
obtaining numerical results. 

Notice that Eq. (82) can now be evaluated without 
performing any quadritures once the appropriate 
solutions to the Rayleigh equation (49) have been 
found. But this equation is still fairly complex and a 
great deal of modeling and/or computation is required 
to evaluate all of its terms. It is therefore worthwhile to 
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introduce some approximate results that can be used to 
get some insight into the import physical effects. 

Some Approximate Results 

It is well known 20 that the Reynolds stresses 7^°* and 

7] y* that appear in the lowest order mean flow 

equations can be chosen so that these equations possess 
similarity solutions of the form 

U = U(ij), R = R(>7) (85) 

where 

/7 = x 2 /x (86) 

and since this causes no particular inconsistency with 
the expansions (33) through (37) and agrees well with 
experimental observation we shall now assume this to 
be the case. 


We also suppose that the observation point lies in the 
plane perpendicular to that of the shear layer, i.e. that 
the azimuthal angle cp = 0. Then we can evaluate all of 
the factors in Eq. (82) by solving the one-dimensional 
Rayleigh equation, i.e. by setting & 3 = 0 (or equivalently 
k = ki) in Eq. (47) which now becomes 


where 


and 




d cl 
d, l {aU-lf 



(O = Xw 


<T = k\! CO 


(aU-\f 


(87) 


(88) 

(89) 


is the reciprocal of the complex phase speed. The 
eigenvalue k in Eq. (76) can now be replaced by 


(7 0 (o)) = — (90) 

co 


i.e. it is a function of the single variable co which 
means that a 0 is now determined by 


cr 0 (q) = — cos# (91) 

c + 

where Q = coa Q is essentially the complex frequency 
eigenvalue of the Rayleigh problem corresponding to 
the specified real phase speed cj cos (9. 

Additional simplification can be achieved by assuming 
that 


( 92 ) 


which implies, among other things, that h' a - 0 .For 
simplicity we also neglect the cross coupling between 

the two components of the Greens function and set c~ 

— c + = constant. It then follows from Eqs. (45), (47), 
(65) , (74) through (77) and (82) that 




co 3 X'( r + 1) 


2 7rc\r 


w_ 

' - 1 o) 

1 J ,co, — cos# 

c + ) 

A 

1 ^ 

co, — cos# 

v. C + J 



® 0 ( <o(l-M 0 cos#))- 


w 0 _{jj',co)a{q\X') 


2 k -(2/)Im©„(0,®) 


Re ® 


sa n 


A 'H 


0 (®(l-f/ c o- 0 (®))) 


where 



[aU-l) 



(j — dcr 


(J 


O ~ 


dc X> 

dn 


(93) 


(94) 

(95) 

(96) 


and W„ denotes the Jacobian with derivatives with 
respect to rj. 

Numerical Results 

To compute the radiated sound field from Eq. (93) it is 
necessary to specify the turbulence spectra ® 0 . Recent 
numerical studies 18 ’ 20 suggest that it should be of the 
Gaussian form 


O 


O 



e -{a>/2o} s f 


(97) 


where w s =X'w s , £ s =l s /X', p s , and U s are 
constants with obvious physical interpretation and a s is 
an empirically determined scale factor. It is also 
necessary to specify the mean velocity profile U ( rj) in 
the shear layer which is reasonably well represented 

by 20 


U = — (l-tanh^) 


(98) 
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The result is, as it should be, independent of the axial s \6 

source location X' within the shear layer when — Re cD 0 (co-U c k[X' , co)^H [-\u\k{X' 

expressed in terms of the scaled frequency co = X'co . v c °° ) 


When Im ©*, =0(1), the 1 Is factor in the exponent of 
Eq. (93) will cause the second term to be negligibly 
small when Im ©*, > 0 and very large (i.e. it will be 
dominant) when Im < O.The s in the denominator 
will cause his term to be large when Im ® rj =0( s )no 
matter what its sign . 


where H(x) denotes the Heaviside function II{x) = 1 ,0 
x 0, is relatively independent of angle and has a much 

narrower width than ® a ( co) due to the relatively narrow 
band of unstable frequencies at X' . 


Concluding Remarks 


Figure 3 is a plot of Im ©„, as a function of 6 for 
various values of the Mach number M and source 
frequency co . The curves are truncated at small 0 
because the second term has no radiation field beyond 
this point, i.e. it becomes evanescent there. These 
results show that Im © w is always positive at 9 = 90° 
but can change sign and become negative when M is 
sufficiently large and 0 is sufficiently small. 

The spectral shape at 90° is therefore produced by the 
first term in Eq. (93) which is primarily determined by 
the factor 



at that angle. Refraction effects, which are accounted 
for by the first factor in this term, will cut off the high 
part of the source spectrum at small angles to the axis. 
But figure 3 shows that the second term can become 
dominant at these angles when M is sufficiently large. 
Its spectral shape, which is primarily determined by 


The Navier-Stokes equations, rewritten in the form of 
the LNS equations with externally applied stress and 
energy flux sources were solved for a two-dimensional 
shear layer by using a vector Greens function approach 
and assuming that the spread rate s of the shear layer is 
small. The relevant Greens function has two 
components -one of which corresponds to the usual 
Lilley equation solutions that appear in the literature 
and the other associated with linear instability waves. 

Numerical results obtained for a simplified source 
model show encouraging qualitative agreement with 
experimental observations. Observed differences in 
spectral shapes at small and large angles to the 
downstream axis that occur at the higher Mach numbers 
can, in part, be attributed to the fact that the first Greens 
function component is dominant at 9 = 90° while the 
second component is dominant at small 6 . 
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